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ABSTRACT 

o , 

We present an interferometric measurement of the Sunyaev-Zel'dovich effect (SZE) 
at 1 cm for the galaxy cluster Abell 2163. We combine this data point with previous 
measurements at 1.1, 1.4, and 2.1 mm from the SuZIE experiment to construct the 
most complete SZE spectrum to date. The intensity in four wavelength bands is fit 
to determine the Compton y-parameter (y ) and the peculiar velocity (v p ) for this 
cluster. Our results are y = 3.56±g;£j:g;^ x 10~ 4 and v p = 410jl^ ±^ km s" 1 
, where we list statistical and systematic uncertainties, respectively, at 68% confidence. 

These results include corrections for contamination by Galactic dust emission. We find 
^r) [ less contamination by dust emission than previously reported. The dust emission is 

distributed over much larger angular scales than the cluster signal and contributes little 
to the measured signal when the details of the SZE observing strategy are taken into 
account. 
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Subject headings: galaxies: clusters: individual (Abell 2163) — Sunyaev-Zel'dovich 
q , Effect — cosmic microwave background — techniques: interferometric 



1. INTRODUCTION 



The thermal Sunyaev-Zel'dovich effect (SZE) is a spectral distortion in the Cosmic Microwave 
Background (CMB) which results from the inverse Compton scattering of low-energy CMB photons 
by thermally distributed hot electrons in massive clusters of galaxies. The thermal SZE manifests 
itself as an increment in the Wien part of the CMB spectrum and a decrement in the Rayleigh- 
Jeans part of the spectrum (Sunyaev & Zel'dovich 1972), with a crossover wavelength of ~ 1.38 mm. 
There is also a kinetic SZE, caused by the bulk radial motion of the cluster with respect to the CMB 
rest frame. This effect shifts the CMB Planck spectrum to a slightly lower (higher) temperature for 
positive (negative) peculiar velocities where positive is defined as motion away from the observer. 
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In this paper, we present a recent 1 cm interferometric detection of the SZE in the cluster 
Abell 2163, obtained with the BIMA and OVRO arrays. We combine this measurement with a 
previous result from the Sunyaev-Zel'dovich Infrared Experiment (SuZIE), a six element bolometer 
array (Holzapfel et al. 1997a, hereafter H1997a) on the 10.4 m Caltech Submillimeter Observatory. 
The SZE in Abell 2163 was also observed by Diabolo, a two element bolometer array mounted on the 
IRAM 30 m telescope (Desert et al. 1998), and their result is included in Table 1. The measurement 
has considerable uncertainty and carries very little weight in the spectral fit relative to the SuZIE 
data. Due to the complication of estimating the effects of systematics on the Diabolo measurement, 
we have chosen to exclude it from our spectral fit. Thus, combining the four measurements from 
SuZIE and OVRO/BIMA, we construct the SZE spectrum and use it to determine the central 
Compton y-parameter and constrain the radial peculiar velocity of the cluster. We also address the 
concern that Galactic dust along the line of sight to the cluster may be a significant contaminant 
in the SuZIE bands (Lamarre et al. 1998). Using data from the Infrared Space Observatory (ISO) 
and the Infrared Astronomical Satellite (IRAS) we determine the level of dust contamination in 
the SuZIE measurements. 

In §2, we describe our observations, data reduction, and dust analysis. In §3 we describe the 
theoretical model for the SZE spectrum and present the results of fitting our data to this model; 
we also consider several systematic uncertainties affecting the fit. All uncertainties are quoted at 
the 68% confidence level unless otherwise specified. 

2. OBSERVATIONS AND DATA REDUCTION 

The 1 cm observations were carried out at the Berkeley Illinois Maryland Association (BIMA) 
millimeter array and the Owens Valley Radio Observatory (OVRO) millimeter array during the 
summers of 1996, 1997, and 1998. The telescopes were outfitted with centimeter-wavelength re- 
ceivers specifically designed to detect the SZE in distant galaxy clusters (Carlstrom et al. 1996). 
The receivers use cooled high electron mobility transistor (HEMT) amplifiers (Pospieszalski et al. 
1995) sensitive to radiation between 0.83 and 1.15 cm with typical receiver temperatures in the 
range 11-20 K. Typical system temperatures (scaled to above the atmosphere) for the Abell 2163 
observations range between 45 and 55 K. 

We observed Abell 2163 at OVRO in 1996 and 1998 for a total of 37 hours with six 10.4 m 
telescopes providing a 4'. 2 FWHM primary beam; we used two 1 GHz wide bands centered at 
1.0 cm and 1.05 cm. This cluster was observed at BIMA in 1996, 1997 and 1998 for 11 hours with 
nine 6.1 m telescopes providing a 6'. 6 FWHM primary beam; we used an 800 MHz bandwidth 
centered at 1.05 cm. At both observatories, we interleaved cluster measurements with observations 
of the radio point source 1549+026 (3.60 ± 0.01 Jy) every 25 minutes to monitor system gains, 
and we used Mars for our absolute flux calibration (Grego et al. 2000). We reduced the data 
using the MIRIAD software package (Sault et al. 1995) at BIMA and the MMA package (Scoville 
et al. 1993) at OVRO. In all observations, we flagged data from baselines where one telescope was 
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shadowed by another, data that were not bracketed in time by observation of a phase calibrator, 
and data showing evidence of poor atmospheric coherence or spurious correlations. Figure 1 shows 
a deconvolved SZE image of Abell 2163 constructed from BIMA data. SZE images are produced 
with the software package DIFMAP (Pearson et al. 1994). In Figure 1 we have applied a Gaussian 
taper with a half power radius of 1.5 kA to the u-v data. A bright point source has been removed 
from the image in Figure 1; see §2.1 for more details. 

An interferometer measures the Fourier transform of the sky brightness, hence we fit a model 
to the data in the Fourier plane where the noise characteristics and spatial filtering of the inter- 
ferometer are well understood. The model is constructed in the image plane, multiplied by the 
primary beam of the telescope, then Fourier transformed to the conjugate plane, known as the u-v 
plane, where it is compared with the data. 

We model the cluster gas as a spherical isothermal-/? model (Cavaliere & Fusco-Femiano 1978), 
where the SZE intensity change is given by 



Here AI U0 is the central SZE decrement /increment at frequency u, 9 C is the angular core radius of 
the cluster, and [3 is a power law index. We adopt the model shape parameters 9 C = 1'.20 ± O'.ll 
and (3 = 0.616 ± 0.031 from an analysis of ROSAT X-ray data (Holzapfel et al. 1997b, hereafter 
H1997b). With these shape parameters we find AI 30 = —0.048 ± 0.006 MJy sr _1 (corresponding to 
a CMB temperature decrement AT CMB = —1780 ± 210 fiK) where A/ 30 is the SZE decrement at 
30 GHz. The statistical uncertainty in the central decrement is determined by varying A/ 30 while 
keeping the shape parameters fixed; all other parameters are free to assume their best-fit values. 
The uncertainties in 6 C and (3 quoted above lead to an additional 0.001 MJy sr _1 uncertainty on 
the 1 cm decrement. 

The 1 cm decrement and combined statistical uncertainty are listed in Table 1, which also 
contains all spectral SZE data for this cluster available to date. The SuZIE AI U0 values are from 
fits to the SuZIE data using the average of the measured spectral responses of the individual 
channels. The AI U0 value from the Diabolo experiment is from Desert et al. (1998). As the table 
indicates, we have made small corrections in the SuZIE bands to account for dust contamination 
along the line of sight to the cluster. These corrections are described in more detail in §2.2. 



Point sources are identified in DIFMAP using a high-resolution map produced with data from 
baselines 20 meters and longer. We find one point source in the Abell 2163 data offset —40" in 
RA and 16" in declination from the SZE pointing center, a position consistent with that of a 
variable radio source with an inverted spectrum detected by the VLA (M. Birkinshaw 2001, private 
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2.1. 



Radio Point Sources 
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communication) at 2 cm and 6 cm. The point source has been removed from the image in Figure 
1. 

The flux and position of the detected point source is jointly fit with the SZE decrement, but 
we must consider the uncertainty in the decrement introduced by point sources in the field with 
flux densities below our detection threshold. Such a point source near the cluster SZE center would 
cause us to underestimate the magnitude of the decrement, whereas a point source in a beam 
sidelobe could lead to an overestimate. To determine the resulting uncertainty we have randomly 
added point sources to the data as follows. Progressing in steps of 1 /xJy, we choose a Poisson 
deviate of the expected number of point sources at each flux density. We add sources to the field 
accordingly out to a radius of 5', well past the half power point of the BIMA beam where faint 
point sources should have a negligible effect on the decrement. After reaching our 3a flux density 
limit, we refit the SZE model without accounting for the added point sources. 

To determine the number of point sources expected at each flux density, we first determine a 
point source number count versus flux density relationship of the form 



from our sample of 41 clusters observed at BIMA, where n is the number density of point sources 
and S is their flux density. We find that the distribution of sources beyond 10" from cluster centers 
agrees well with the expected distribution from a larger sample of point sources at 31 GHz (C. Pryke 
2002, private communication). Inside 10", however, we detect nearly 30 times more point sources 
than expected due to the strong contribution of cluster cD galaxies. We therefore exclude the inner 
10" of the map when selecting point sources to construct dn/dS and address the issue of central 
point sources in more detail below. We choose 240" as an outer limit in selecting point sources 
as it is a good trade-off between encompassing a large number of detected sources and allowing 
a low beam-corrected flux density threshold. We use only point sources with beam-corrected flux 
densities greater than 2 mJy, since a 2 mJy point source would have a beam-attenuated flux density 
of 0.75 mJy at 240" and 0.75 mJy is a good fiducial 3a threshold for all 41 maps. Lastly, we assume 
the power law index determined from the Pryke sample, a = 2.4, in Equation (2). Solving for the 
normalization, we find 



We then randomly distribute point sources according to this empirical relationship, starting from 
a flux density of 10 //Jy (tests of the Monte-Carlo showed that starting at a lower flux density has 
a negligible effect on the decrement) up to a flux density of 300 fiJy, corresponding to the 3a limit 
of the more sensitive OVRO Abell 2163 data. The resulting distribution of best-fit decrements is 
nearly Gaussian and centered at the original decrement with a standard deviation of 0.002 MJy sr _1 . 
We adopt this as the statistical uncertainty due to randomly distributed undetected point sources, 
and it has been combined in quadrature with the other statistical uncertainties in Table 1. 



dn/dS oc S 01 mJy 1 arcmin 
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We attempt to estimate the uncertainty from undetected point sources a second way, using 
information on the distribution of point sources in the field from observations at longer wavelengths. 
Sources with flux densities greater than ~ 2 mJy at 21 cm appear in the NVSS catalog (Condon 
et al. 1998). The catalog contains four point sources within 240" of Abell 2163 that we do not detect 
in our 1 cm observations. Searching the 100 /jJy rms OVRO map at the positions of these sources, 
we find no evidence of emission from three of the four; their positions correspond to signals at the 
0.2<7, — 0.4<7, and —la levels. The fourth, however, displaced 30" in RA and —127" in declination 
from the map center, corresponds to a 2<7 signal in the OVRO map. We do not marginalize over 
the fluxes of all four sources because of the loss of degrees of freedom in the fit, but marginalizing 
over the 2a source we find no change in either the decrement or the uncertainty in the decrement, 
which is dominated by the noise in the map. 

To account for the fact that on-center point sources are statistically more common we place a 
point source model at the estimated position of Abell 2163's bright central galaxy (RA 16^15 m 49 s .O, 
Dec — 6°8'42".7 from Digitized Sky Survey) and marginalize over its flux. The new best-fit decre- 
ment is —0.049 MJy sr^ 1 , an increase in magnitude of just 0.001 MJy sr _1 from the original decre- 
ment. The new decrement falls easily within the 68% uncertainty of the original; we revisit this 
result when we discuss systematic uncertainties (§3.1). 

2.2. Analysis of Dust Emission 

Lamarre et al. (1998) (hereafter L1998) determined that IR cirrus flux is a non-negligible con- 
taminant in the three SuZIE bands by including submillimeter data from the PRONAOS balloon 
experiment (Serra 1998) in their fit to the SuZIE and Diabolo data,. They accounted for dust 
by fitting to a combined SZE-spectral model and modified blackbody model, and found that re- 
moving the dust signal changed AI UQ significantly. Specifically, they found that the 1.1 mm signal 
shrinks by 0.051 MJy sr™ 1 to 0.236 MJy sr _1 , the 1.4 mm signal changes by 0.018 MJy sr^ 1 to 
-0.124 MJy sr -1 , and the 2.1 mm band changes by 0.004 MJy sr -1 to -0.387 MJy sr -1 . L1998 
assumed in their analysis that the SuZIE scan pattern is identical to the PRONAOS chopping 
scheme. However, PRONAOS performed a 6'. 2 chop at constant elevation (see L1998 for details) 
while SuZIE uses a drift scan at constant declination (H1997b). Thus, the experimental filters are 
different and will not necessarily detect the same dust signal. 

To better estimate the level of dust contamination in the high frequency SuZIE channels, we 
have applied the SuZIE observing scheme to extrapolated dust maps at the SuZIE wavelengths. 
The expected dust contamination is listed in Table 2 and the dust-corrected SuZIE measurements 
are given in Table 1. Here we describe the procedure. 

For clarity, we summarize the map construction procedures as applied to the 1.1 mm band; the 
1.4 mm and 2.1 mm procedures are identical. The 1.1 mm dust map is created by extrapolating 
from three existing dust maps centered on Abell 2163; ISO maps at 90 /xm and 180 [im, and an 
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IRAS map at 100 /im. We extract a 15' x 1'.75 central region from each dust map, consistent with 
the area of the SuZIE scans. We then smooth each map to a resolution consistent with the SuZIE 
beam size (1'.75), taking into account the intrinsic resolution of the IRAS and ISO experiments. 
We use a fiducial pixel to pixel uncertainty of 0.2 MJy sr~ x (cf. Beichman et al. 1988) for the IRAS 
map; for the ISO maps, we adopt the values estimated in L1998 of 0.14 MJy sr _1 and 0.18 MJy sr _1 
at 90 and 180 /xm, respectively. For each ISO map, we then bin the 0'.25 pixels into 0'.75 pixels, 
consistent with binning performed in the SuZIE analysis. For each new 0'.75 pixel, we take the 
average of the original three 0'.25 pixels as the intensity. The IRAS pixels are roughly twice the size 
of the binned ISO pixels, so we linearly interpolate to estimate the intensity of interleaved pixels, 
which we artificially add to the IRAS map. We then fit a modified blackbody to these dust maps 
at each pixel position; since we have three dust maps, each fit is to three data points, or one "pixel 
triplet". This allows us to extrapolate the intensities out to 1.1 mm. We fit a modified blackbody 
of the form 

F - * e hv/kT d _ V ( 4 ) 

where 7 is the spectral index, Td is the dust temperature, and we normalize the model with the 
dust flux at 180 //m, ^(180). Following L1998, we fix the spectral index at 7 = 2. Upon fitting 
each dust pixel triplet to this model, we find i^(180) in the range 25.1 to 26.1 MJy sr _1 and in 
the range 21.3 to 21.9 K over all triplets. We then use each triplet's best fit i^(180) and Td with 
Gaussian error propagation to extrapolate the intensity at that position and create our 1.1 mm 
dust map. 

After repeating this process at 1.4 and 2.1 mm, we perform the simulated scan following the 
detailed discussion of SuZIE's instrumental specifications and observing scheme in H1997b. In our 
simulation we scan across the 0'.75 pixels of the dust maps in the direction of increasing right 
ascension (RA). For each pixel along the scan, we record the difference between its intensity and 
that of the pixel advanced 4'. 6 in RA; pixel-to-pixel uncertainties are combined in quadrature. This 
procedure is equivalent to SuZIE allowing the sky to drift overhead with the earth's rotation and 
recording the difference in signal between two bolometer array elements separated by 4'. 6 on the 
sky. We do not treat SuZIE's 2'. 3 "triple beam chop" here as it has little sensitivity to linear 
changes in brightness across the sky and therefore is less sensitive to dust; the results presented 
here serve as an upper limit to the dust contribution (see H1997b for more details). Our simulation 
thus creates "differenced" maps of the dust brightness. 

We perform one-dimensional x 2 fits in the parameter A/j, to the dust maps, fitting to a SuZIE- 
differenced (3 model with the pixel-to-pixel dust noise as our uncertainty. We use a Monte Carlo 
procedure and successively fit to 10 6 different realizations of the dust; in a single realization, each 
pixel has added to it some Gaussian-distributed amount of the pixel-to-pixel noise. The results 
of the Monte Carlo are highly Gaussian as expected; we thus take the most frequently returned 
AI U0 as our best fit and those that enclose 68.3% of the returned A/^ 's as our la uncertainties. 
The results are shown in Table 2 and reflected in the last column of Table 1. We find a dust 
contamination level roughly seven times smaller than that reported in L1998, and small compared 
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to the statistical uncertainties of the SZE measurements. 



3. SZE SPECTRAL MODEL AND RESULTS 

We have fit the dust-corrected SZE spectral data in Table 1 to a model which consists of both 
the thermal and kinetic SZE components. The thermal component, A/ T , can be written as 

AI T = I yf(x)(l + 8 T ) (5) 

where x = hv /kT CMB , T CMB = 2.728 K is the CMB temperature (Fixsen et al. 1996), 5 T is a 
relativistic correction to the thermal effect good to fifth order in kT e /m e c 2 (Itoh et al. 1998), 
J = 2{kT CMB ) 3 /(hc) 2 , and 



x 4 e x \x{e x + 1) 



(e x - 1) 



1 



(6) 



The quantities h, k, and c are the Planck constant, Boltzmann constant, and speed of light, 
respectively. The quantity y is the Compton y-parameter, and is given by 



kT e 
nipC 2 



n e a T dl, (7) 



where T e is the electron temperature of the intracluster medium (ICM), m e is the electron mass, 
n e is the electron density, and ax is the Thomson scattering cross section. The Compton y- 
parameter is proportional to the pressure integrated along the line of sight (dl) through the cluster; 
we parameterize it in terms of a central value, which we denote y . 

The kinetic component, AI K , depends on both y and the radial peculiar velocity of the cluster, 
v p , through the equation 

AI K = -{m e c 2 /kT e )I y^g(x)(l + S K ) (8) 

where g(x) = x A e x /{e x — l) 2 and S K is a relativistic correction to the kinetic effect good to third 
order in kT e /m e c 2 (Nozawa et al. 1998). Corrections of order v p 2 /c 2 to the thermal and kinetic 
SZE are neglected since v p <^ c. The total intensity is then the sum of the thermal and kinetic 
components, 

A/„ = A/ T + AJ K . (9) 

We perform a least squares fit to this model in the parameters y and v p assuming an ICM tem- 
perature T e of 12.4+ 2 ! keV from H1997b. 

The resulting v p and y from the least squares fit to the four data points are y = 3.56 ± 0.41 x 
10 -4 and v p = 410^g5Q° km s _1 with x 2 = 0.8 for two degrees of freedom. We compare this with 
the results of our fit to the SuZIE data only, y = 3.72 ± 0.51 x 10~ 4 and v p = 2901^° km s" 1 . 
The addition of the 1 cm point improves the y constraints by 18%. The constraints on v p hardly 
change when the 1 cm point is included since constraints on peculiar velocity are dominated by the 
SuZIE data points which straddle the thermal SZE null. 
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Figure 2 shows the SZE spectrum generated from the fit to all four frequency bands. There is 
excellent consistency amongst all of the measurements. The data were obtained from independent 
experiments, each with its own telescope, detector, and observing strategy, operating at four dif- 
ferent wavelengths. The consistency amongst the data attests to the reliability and accuracy with 
which the SZE is now being measured. 

3.1. Systematic Uncertainties 

In this section we address several additional sources of uncertainty in our spectral fits; the 
results of this analysis are summarized in Table 3. 

We first consider the systematic uncertainty associated with the dust emission calculated in 
§2.2. When we alter the dust intensity at the 68% level, the spectral fit results change by ^.q'qi x 10~ 4 
for y and +g km s _1 for v p . The dust uncertainties are sufficiently small as to have very little 
impact on the spectral fit results. We find a more significant change when we repeat the dust 
analysis with a smaller spectral index, 7. Finkbeiner et al. (1999) recently suggested based on fits 
to FIRAS data, that while 7 = 2 is appropriate for wavelengths < 600 /im, the spectral index at 
longer wavelengths is closer to 1.5. Setting 7 = 1.5 for the entire dust spectrum we find that y 
increases by 0.02 x 10~ 4 and v p decreases by 70 km s _1 . We combine these in quadrature with the 
changes that resulted from scaling the dust at the 68% level and adopt the total as the systematic 
uncertainty due to dust contamination. 

Temperature fluctuations in the CMB are a potential source of confusion for SZE measure- 
ments. Measurements of the kinetic SZE are particularly susceptible since the spectral dependence 
of the kinetic SZE and CMB fluctuations are identical, absent relativistic corrections to the SZE 
at the few percent level. Recent BIMA observations at 28.5 GHz have provided good constraints 
on CMB anisotropy at the angular scales of our measurement (Dawson et al. 2001; Holzapfel et al. 
2000); they find AI V < 4.8 x 10~ 4 MJy sr" 1 at 95% confidence, which is just 1% of A/ 30 . This 
uncertainty is negligible given the precision of our measurement. The level of CMB contamination 
may be larger in the SuZIE bands, however, particularly near the null in the thermal SZE. 

To test the level of CMB contamination in the SuZIE bands we first use the software package 
CMBfast (Zaldarriaga Sz Seljak 2000) to generate a CMB power spectrum up to a multipole £ = 
2500. We assume a universe consisting of cold dark matter with a cosmological constant (ACDM) 
and use conventional parameter values Sldm, r c , n s , h) = (0.04, 0.26, 0.70, 0, 0.95, 0.72) where 
Clb is the fractional baryon density of the universe, fl c dm is the fractional cold dark matter density, 

is the fractional density of vacuum energy, r s is the optical depth to reionization, n s is the 
spectral index of the power spectrum, and Hq = 100/i km s _1 Mpc -1 is the Hubble constant. The 
values chosen are good estimates given the range of values currently constrained by CMB anisotropy 
measurements (Pryke et al. 2001; Netterfield et al. 2001; Stompor et al. 2001); h comes from the 
Hubble Space Telescope Hubble Constant key project (Freedman et al. 2001). To account for sec- 
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ondary anisotropies we conservatively assume a flat power spectrum for £ > 2500, corresponding 
to arcminute angular scales. The flat band power we use corresponds to AT CMB /T CMB ~5x 10~ 6 , 
which is consistent with expectations of secondary anisotropies (Hu & Dodelson 2001) and the 
BIMA limit of Dawson et al. (2001). We then use this power spectrum to produce mock CMB 
fields at SuZIE's resolution. After simulating the SuZIE scan on these maps, we fit them to an 
isothermal /3-model and use a Monte Carlo procedure similar to that used with the dust to deter- 
mine the contribution of CMB to the total signal in each SuZIE band. Upon fitting to 10 4 different 
maps for each band, we find an rms of ±0.020 MJy sr _1 at 1.1 mm, ±0.022 MJy sr -1 at 1.4 mm, 
and ±0.017 MJy sr _1 at 2.1 mm. Further tests using only primary anisotropy {I < 2500) or only 
secondary anisotropy (£ > 2500) indicate that primary anisotropy is responsible for roughly 80% 
of the total CMB contribution in all three bands. Including the rms total anisotropy contributions 
in the SuZIE signals listed in Table 1, we redo the spectral fit and find changes of ±0.05 in y 
and ^^qq km s _1 in v p . The peculiar velocity suffers more from confusion with CMB anisotropy 
than does the Compton y-parameter, as anticipated. We adopt these results as the systematic 
uncertainties due to CMB contamination. 

We also consider confusion with extended emission from the radio halo surrounding Abell 2163, 
which is the brightest such halo yet discovered (Feretti et al. 2001). To test its effect on the 1 cm 
decrement, we combine the OVRO data with the 1.4 cm NVSS map of the cluster region, which 
shows a large region of extended emission surrounding the cluster. We first scale the NVSS map 
to account for the unknown spectrum of the radio halo; we make this scale factor a free parameter 
and redo the fit described in §2. The decrement is increased by —0.003 MJy sr -1 when we include 
this extra parameter, and the scaling is consistent with a spectral index a ~ 1.2 where the intensity 
has a v~ a dependence. We note, however, that our fit is consistent with a = oo (i.e., no radio halo) 
within 68% confidence. Furthermore, the spectral index may not be constant across the halo, so 
the best-fit a must be taken as a crude estimate. H1997b also address radio halo confusion, and 
estimate a 2% effect at 2.1 mm. This assumes a spectral index of ~ 1.5, however; using our best-fit 
value a = 1.2, the effect increases to ~ 5%. When we alter the 1 cm and 2.1 mm points accordingly 
and redo the spectral fit, we find that y increases by 0.17 x 10 -4 to 3.73 x 10 -4 , and v p increases 
by 50 km s _1 to 460 km s _1 . We adopt these as estimates of the systematic uncertainties in y and 
v p due to radio halo confusion. 

We use a similar treatment to determine the systematic uncertainty due to point source con- 
tamination. In §2 we accounted for the possibility of an undetected on-center point source in 
the field at 1 cm, and found that it could increase the magnitude of the measured decrement by 
0.001 MJy sr" 1 . Such a faint source would not affect the SuZIE data unless it had a strongly in- 
verted spectrum; we assume that it does not. The brighter off-center point source detected at 1 cm 
could affect the SuZIE data if its spectrum were inverted at millimeter wavelengths. Although its 
flux rises from 1 to 3 mJy between 6 cm and 2 cm (M. Birkinshaw 2001, private communication), 
we found at the time of our observations that the flux of the source at 1 cm varied between 1.1 and 
1.4 mJy. This corresponds to a spectral index a between 1.1 and 1.5. Extrapolating with such a 
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spectral index to 2.1 mm, we find a flux of less than 0.2 mJy which would have a negligible effect 
on the decrement. However, because of its variability and unusual spectrum, we conservatively 
assume a flat spectral index (a = 0) for A < 1 cm. Changing the 1 cm decrement to account for a 
faint central point source and altering the SuZIE 2.1 mm decrement to account for a 1.5 mJy point 
source, we find an increase of 0.07 x 10~ 4 in y and 50 km s _1 in v p . We include these in Table 3 
as our systematic uncertainty due to point source contamination. 

We test the importance of our assumed 12.4 keV electron temperature to y and v p by per- 
forming the spectral fit over a large range in T e . An increase in T e trades off with a decrease in 
the cluster optical depth r = J n e aTdl, as can be seen in equation (7); thus so long as r remains 
a free parameter, Compton y is only weakly sensitive to temperature changes. In the la range of 
T e (10.5 keV < T c < 15.2 keV), y varies by j^'o^ x 10~ 4 . The peculiar velocity also has a weak 
temperature dependence; it rises with temperature at low T c (< 5 keV) as r decreases (cf. equation 
8), then falls at high T e (> 15 keV) due to relativistic corrections. The la range of T e corresponds 
to a broad plateau in the T e -v p curve, and we find v p varying by only +20 km s _1 and —40 km s -1 . 

We have assumed an isothermal temperature distribution for the cluster gas, consistent with 
the H1997a and b treatments of Abell 2163. Although the assumed cluster temperature has little 
impact on y and v p , uncertainties in the thermal structure of the cluster may have a larger effect. In 
a recent study of Abell 2163 using the BeppoSAX X-ray satellite, Irwin & Bregman (2000) find the 
cluster temperature profile consistent with isothermality at the 2a confidence level out to a radius 
of 9'; this is consistent with the result of White (2000), who finds a temperature profile consistent 
with isothermality at the la confidence level out to 10'. Chandra or XMM data are necessary 
to place stronger constraints on the thermal profile of Abell 2163, and we postpone quantifying a 
systematic uncertainty until such data become available. 

Finally, we consider the effect of the assumed ±8% absolute calibration uncertainty in each 
SuZIE band and the ±3% absolute calibration uncertainty in our interferometric 1 cm measure- 
ments. These combine to change y by as much as ±0.17 x 10~ 4 and v p by up to ±170 km s _1 . 

4. CONCLUSIONS 

We have used interferometric data obtained at the OVRO and BIMA observatories to deter- 
mine the SZE intensity decrement in the cluster Abell 2163 at a wavelength of 1 cm. We find 
A/ 30 = -0.048 ± 0.006 MJy sr" 1 by modeling the ICM as an isothermal spherical /3-model. We 
have thus expanded the spectral coverage of this cluster to include the Rayleigh-Jeans part of 
the CMB blackbody spectrum and we present the most complete SZE spectrum to date. Upon 
fitting the OVRO/BIMA and SuZIE data for Abell 2163 to an SZE spectral model, we find 
y = 3.56toi\toil x 10~ 4 and v p = 4W±Hl°±f^ km s _1 where we list statistical uncertainty 
followed by systematic uncertainty at 68% confidence. The systematic uncertainties in Table 3 
have been combined in quadrature to determine the totals reported above. We have shown that 



- 11 - 



adding the 1 cm measurement improves constraints on the Compton y-parameter, as expected, but 
has less effect on the peculiar velocity which is best constrained by measurements near the 218 GHz 
thermal null. 

The spectral fit includes corrections for dust emission in the SuZIE bands; we find that the 
contamination level depends strongly on the observing scheme of the instrument due to the distri- 
bution of dust on the sky. The dust contamination level for a cluster will obviously depend on its 
location on the sky; the line of sight to Abell 2163 passes near the Galactic plane, so it is more 
susceptible to dust contamination than clusters at higher galactic latitudes. 

Finally, we emphasize that the Abell 2163 spectrum shown in Figure 2 consists of data from 
two very different experiments; an array of six bolometers performing drift scans (SuZIE) and two 
interferometers (OVRO and BIMA). The remarkable consistency that we find amongst the two sets 
of measurements is further evidence that SZE observations are now a reliable probe of structure in 
the universe. 
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Right Ascension (J2000) 

Fig. 1. — Deconvolved 1 cm SZE image for Abell 2163, constructed from BIMA data. The image 
was made by first applying a Gaussian taper to the u-v data with a half power radius of 1.5 kA; 
the resulting beam FWHM (110" x 90" at position angle 30°) is shown by the white ellipse in the 
lower left. The rms is 300 /xJy beam -1 , corresponding to a Rayleigh Jeans brightness sensitivity of 
50 fiK. Contours are integer multiples of ±1.5cr with negative contours shown as solid lines. 
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Fig. 2. — SZE spectrum of Abell 2163 (points) and best fit model (lines). The dashed line is the 
thermal spectrum, dot-dashed the kinetic spectrum, and solid the sum of the two; SuZIE data 
points appear as triangles and OVRO & BIMA as a circle. The insert shows the 1 cm point with 
its error bar. 



Table 1. Measurements of SZE in Abell 2163 



A 


Instrument 


Measured AI V0 


Dust-corrected AI U0 


(mm) 




(MJy sr" 1 ) 


(MJy sr" 1 ) 


10.0 


OVRO & BIMA 


-0.048 ± 0.006 




2.1 


Diabolo 


-0.545 ±0.22 




2.1 


SuZIE 


-0.381 ± 0.037 


-0.380 ± 0.037 


1.4 


SuZIE 


-0.106 ±0.077 


-0.103 ±0.077 


1.1 


SuZIE 


0.287 ±0.105 


0.295 ±0.105 
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Table 2. Strength of Dust Signals 



Map 

(H 


(MJy sr" 1 ) 


1100 
1400 
2100 


-7.6±i;4 x 10~ 3 
-3.51°;^ x 10" 3 
-0.7j£J x 10~ 3 



Table 3. Systematic Uncertainties in y and v p 



Source 


Effect on y 


Effect on v p 




(xlO- 


- 4 ) 


(km s" 




Dust Contamination 


(+0.02,- 


-0.01) 


(+30, 


-70) 


CMB Anisotropy 


(+0.05,- 


-0.05) 


(+420, - 


-400) 


Radio Halo 


(+0.17,-0) 


(+50,0) 


Point Sources 


(+0.07, -0) 


(+50 


,-o) 


X-ray Temperature 


(+0.09,- 


-0.07) 


(+20, 


-40) 


Absolute Calibration 


(+0.17,- 


-0.17) 


(+170,- 


-170) 


TotaP 


(+0.27,- 


-0.19) 


(+460, - 


-440) 



a combined in quadrature. 



